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1. Introduction 

Spin symmetry lias played an important role in the study of heavy hadrons during the 
past several years. Many of its consequences for mesons and baryons with quark content 
Qq and Qqq have been investigated within the context of the Heavy Quark Effective 
Theory (HQET) [|I|. The spin and flavor symmetries of QCD that become exact in the 
Mq -^ oo limit are transparent in this effective theory framework. But even without the 
HQET apparatus, several aspects of charm and bottom hadron phenomenology may be 
understood from symmetry considerations alone. For example, the rates for strong and 
electromagnetic transitions between states belonging to separate spin doublets are simply 
related by group theory factors. Moreover, the sums of partial widths for each member of 
a doublet to decay to states in another doublet are equal. These predictions that follow 
from spin symmetry in the infinite mass limit appear to be reasonably well satisfied for 
finite mass hadrons in the charm and bottom sectors 0. 

In this letter, we consider the implications of spin symmetry for heavy quarkonia. 
Since the spins of both the quark and antiquark inside a QQ bound state are separately 
conserved as Mq — > oo, transitions between members of different quarkonia multiplets 
within the same flavor sector are more tightly constrained than their heavy-light meson 
analogues. Ratios of decay rates can be predicted and checked against the large body of 
charmonia and bottomonia data collected over the past two decades. Spin symmetry also 
constrains QQ bound state formation and ties in closely with recent perturbative QCD 
calculations of parton fragmentation to heavy quarkonia |0-|^ . As we shall see, the union 
of these ideas yields several interesting predictions which can be experimentally tested at 
hadron colliders. 

Our paper is organized as follows. We flrst review the general relation that spin sym- 
metry imposes upon heavy hadron transitions in section 2. We then test this relation for 
quarkonia by comparing its predictions for QQ electromagnetic transitions with experi- 
mental data. In section 3, we use the spin symmetry relation along with nonrelativistic 
QCD power counting arguments to determine the leading order alignment of P and D-wave 
ort ho quarkonia produced via gluon fragmentation. Finally, we discuss the implications of 
spin symmetry for a fragmentation mechanism which may resolve the current large dis- 
crepancy between theory and experiment in prompt tj)' production at the Tevatron. 
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2. Radiative Quarkonia Transitions 

It is useful to recall the basic symmetry relation which constrains transitions between 
any initial and final hadronic states H and H' that contain heavy constituents . For Qq 
mesons and Qqq baryons, we let j)^ , j^ and J^'' denote the angular momentum of H^''^s 
heavy, light and total degrees of freedom. As the mass of the heavy constituents inside H^'' 
tends to infinity, the separate angular momenta of the light and heavy degrees of freedom 
become good quantum numbers. A similar phenomenon occurs for QQ quarkonia. But 
in this case, the quantum number jf^ represents just the total spin of the QQ pair, while 
j) stands for the bound state's remaining angular momentum which includes the orbital 
contribution from the heavy quark and antiquark. 

In the infinite mass limit, the amplitude for a general H —^ H' process that is mediated 
by an interaction Hamiltonian Ti^ and results in the emission or absorption of some light 
quanta can be decomposed as 

lA = {J'J'^IHjIJJ^) 

= ^{J'J'z iJ'eJiz ; J'hJ'hz) iJiJ'ez ; J'hJ'hz I T^i \Jdez ; jhjhz ) {jdiz ; jhjhz \JJz) (2.1) 
= Xl^^'^^ iJefizdhJ'hz) ifefiz \^i \JeJ£z){Jdez;JhJhz\JJz)Sjj'Jj^j'^^ 

where repeated angular momentum labels are summed. If the reaction proceeds through 
the L-th multipole of Hi, the amplitude becomes proportional to the reduced matrix 
element (j^||L||j£) and may be conveniently rewritten in terms of a 6j-symbol [0,0]: 



i^ = \Y^ ( J' J'z iJ'eJiz ; JhJhz ){LL^; jgj'f,^ \jeji^ ) {jejg^ ; jhjhz \JJz) {3'(\\L\ \ji) 

= (_l)^+^;+^'^+V(2j. + 1)(2J' + 1) {^^^ ^j ^^,}(L(J,-j;);J'j;|JJ,)0-;||L||j,). 

(2.2) 
The consequences of spin symmetry follow from this last formula. EJ 

We will apply eqn. (|2.2|) to heavy quarkonia transitions. In order to test its validity, 
we first investigate charmonia and bottomonia radiative decays. We start with the electric 



The only 6j-symbols which we will need for apphcations of (2^) in this letter are 



L L 1 ^-i)L+3h+J ( L L 1 ^-i)L+jh+J' 

~ and 



jh J jh j ^{2L + l){2jh + l) \jh jh J'] ^{2L + l){2jh 



dipole process ^'^Pj -^ ^'^Sj' + 7. In the initial state, the heavy spin angular momentum 
jh = or jh = 1 can be combined with the P-wave orbital angular momentum ji = 1 
to form a /iq or Xqj quarkonium. After the L = 1 radiative transition, the final state 
contains an S-wave tjq or i/jq with j'g = 0. Eqn. ( p.2| ) implies that the decay rate Fi^'^Pj —>■ 
^'^Sj' +7) oc E^, averaged and summed over initial and final polarizations, is independent 
of J and J' in the infinite mass limit. For finite mass quarkonia, the rate indirectly 
depends upon these angular momentum quantum numbers through formally subleading but 
phenomenologically important splittings between members of the S and P-wave multiplets. 
Taking these splittings into account, we find the following radiative partial width ratios in 
the cc sector: 

r(xcO ^ J/ip + 7) : r(xci ^ J/ip + 7) : r(xc2 -> J/ip + 7) : r(/i, ^Vc + l) 

= 0.095 : 0.20 : 0.27 : 0.44 (Theory) (2.3a) 

= 0.092 ± 0.041 : 0.24 ± 0.04 : 0.27 ± 0.03 : unmeasured. (Experiment) (2.36) 

Analogous spin symmetry predictions in the bb sector cannot be checked against data, for 
absolute radiative partial widths have not yet been measured. However, we can compare 
theoretical and experimental values for the branching fraction ratio 



^ r(x..(2P)^T(l^)+7) ^ E,i2P^lS) 
' T{xBj{2P)^T{2S)+^) [E^{2P^2S)i 

as a function of J: 



(2.4) 



Ro:Ri:R2 = 3.22 : 2.55 : 2.28 (Theory) (2.5a) 

= 5.11 ± 4.13 : 2.47 ± 0.60 : 2.28 ± 0.47. (Experiment) (2.56) 

We do not include the corresponding prediction for r{hb{2P) -^ r]i,{2S) + '^)/r{hh{2P) -^ 
r]b{lS) +7) in eqn. (|2.5a|) since the parabottomonium phase space factors are uncertain as 
the rjbilS), r]i,{2S) and hh{2P) masses have not been experimentally measured. 

We can further test spin symmetry for heavy quarkonia by considering ^'^Sj -^ ^'^Pj' + 
7 electric dipole decays. Inserting ji = 0, jg = 1 and L = 1 into amplitude (|2.2|), we find 



We use the experimentally measured paracharmonia masses Mr,^ = 2979 MeV B, Mj 



he 



3526 MeV [|lO|,|l| and M^,^ = 3592 MeV [H to predict the V{hc ^ ??c + 7) and Vii]'^ ^ /ic + 7) 
entries in eqns. ( p. 3a ) and (2.6a). 



T{^Si -^ ^Py + 7) oc ((2J' + l)/9) E^ and T{^So -^ ^Pi + 7) oc E^. These partial width 
results yield the following charmonium and bottomonium sector ratios: 

r(^' ^ xco + 7) : r(^' -> xci + 7) : r(^' -. Xc2 + 1) ■ vW, ^ /ic + 7) 

= 9.3 : 7.9 : 5.4 : 1.3 (Theory) (2.6a) 

= 9.3 ± 0.8 : 8.7 ± 0.8 : 7.8 ± 0.8 : unmeasured (Experiment) (2.66) 

r(T(2^) -^ xm{lP) + 7) : r(T(25) ^ Xhi{lP) + 7) : r(T(25) ^ xUlP) + 7) 
= 4.29 : 6.70 : 6.59 (Theory) (2.7a) 

= 4.3 ± 1.0 : 6.7 ±0.9: 6.6 ±0.9 (Experiment) (2.76) 

r(T(3^) -^ X6o(2P) +7) : r(T(2^) ^ xw(2P) +7) : r(T(25) ^ X62(2P) ±7) 
= 7.0 : 11.3 : 12.3 (Theory) (2.8a) 

= 5.4 ± 0.6 : 11.3 ±0.6: 11.4 ±0.8. (Experiment) (2.86) 

Again we do not include parabottomonium partial widths in eqns. ( p.7a|) and ( p.8a|) since 
their symmetry breaking phase space factors are uncertain. 

Comparing the theoretical and experimental numbers in eqns. (|2.3|) - (|2.8| ), we see 
that the predictions of spin symmetry agree quite well with the quarkonia data. The 
discrepancies in a given 66 mode are smaller than those in its cc analogue, as one would 
expect. Spin symmetry also works best for the lowest lying QQ states where the ratio of 
the heavy quarks' kinetic energy to their rest mass is minimal. But the general success of 
these radiative transition findings bolsters one's confidence that results for other processes 
which follow from the amplitude expression in eqn. ( |2.2| ) will also work well. We therefore 
use it to investigate quarkonia production in the following section. 

3. Orthoquarkonia Alignment 

Quarkonia physics involves a number of different scales separated by the velocity v of 
the heavy quarks inside QQ bound states. These scales are set by the quarks' mass Mg, 
typical momentum MqV and kinetic energy Mqv"^. Within the past few years, an effective 
field theory formalism called Nonrelativistic Quantum Chromodynamics (NRQCD) has 
been established to keep track of this quarkonia scale hierarchy |T^ . The effective theory 



is based upon a double power series expansion in the heavy quark velocity and strong 
interaction fine structure constant. Its power counting rules allow one to methodically 
determine the most important contributions to a particular quarkonium process |T^. We 
will utilize these NRQCD rules in conjunction with spin symmetry arguments to inves- 
tigate quarkonia production via gluon fragmentation. As we shall see, these tools yield 
phenomeno logically interesting predictions for ort ho quarkonia alignment. 

We first consider g -^ Xqj fragmentation. The Feynman diagrams that mediate this 
process at lowest order in the the velocity expansion are illustrated in figs, la and lb. The 
graphs contribute to Xqj production through scattering collisions like gg -^ gg* -^ ggXgj- 
Such fragmentation reactions dominate over lower order parton fusion processes when the 
lab frame energy qq of the fragmenting gluon g* is large, but its squared four-momentum 
q^ is close to the bound state's squared mass M^ ~ (2Mq)^. In this case, g* is transverse 
up to corrections of order q^ /q^- Consequently, the fragmenting gluon is essentially real 
so long as its energy q^ is significantly greater than M^^j. 

The QQ states that emerge from the shaded ovals in fig. |1| reside within the first two 
Fock components of the Xqj wavefunction 

\Xcj) = 0{l)\QQ['pi''^]) + 0{v)\QQ['s['^]g) + ---. (3.1) 

The angular momentum quantum numbers of the QQ pairs are indicated in spectroscopic 
notation inside the square brackets, and their color configurations are labeled by singlet or 
octet superscripts. In fig. |l|a, the incoming hard gluon with four-momentum q fragments 
into QQ[^Pj ^] plus an outgoing hard gluon. The rate for this process is proportional 
to the square of the first derivative of the P-wave bound state's wavefunction evaluated 
at the origin | -Ri(0)| ^ ~ v^. In fig. |l]b, the incoming gluon turns into QQ[^Si ] which 
appears in the next-to- leading Fock component in ( |3.1|) . The formation rate for the colored 
S-wave pair is proportional to the square of its wavefunction | -RslO)] ^ ~ v^. The ^S^ 
state subsequently converts to a Xqj through a soft gluon coupling which is depicted as 
a "black box" in fig. |lp3. This "black box" may be regarded as mediating the dominant 
chromoelectric dipole transition '^S[ -^ "^PJ + g. Although the resulting P-wave state 
has order unity overlap with the full Xqj wavefunction, the soft gluon emission costs one 
power of V in the amplitude and two in the rate. The total transition in fig. |l|b therefore 
proceeds at 0{v^). Alternatively, the "black box" may be viewed as combining the ^Sl 
QQ pair with a soft gluon to form the complete | QQ [^S^ ] g) Fock component in eqn. 



which has 0{v) overlap with | Xqj)- This interpretation again leads one to conclude that 
the transition in fig. |l|b takes place at 0{v^) 0. 

Although the Feynman diagrams in fig. ^ are the same order in the NRQCD velocity 
expansion, their dependence upon the short distance fine structure constant is different. 
The graphs in fig. |l|a involve two hard gluons and are suppressed by as{2MQ)/7T relative to 
the graph in fig. |l|b which has only one. This perturbative QCD suppression is partly offset 
by the numerical values Hi ^ 15 MeV and H^ ~ 3 MeV for the nonperturbative matrix 
elements associated with the color-singlet and octet graphs [§]. However, prompt Xqj 
production at high transverse momenta is known to be dominated by gluon fragmentation 
to QQ['^S[ ] [|T5H T7|. Therefore, we will simply neglect the subleading contributions to 
g ^ Xqj which come from the diagrams in fig. |l|a. El 

The intermediate ^S"} state in fig. |T]b has the same quantum numbers as the incoming 
hard gluon. In particular, it is transversely aligned at high energies. Spin symmetry fixes 
the degree to which the final Xqj inherits the color-octet state's alignment through the 
soft gluon processes "^S"} ^ ^ ^Pj + g or ^Sl + g ^>- '^PJ . The populations of all the 
individual Xqj helicity components may be determined by setting J£ = 0, j^ = 1 and L = 1 
in the general amplitude expression ( |2.2|) . After a straightforward computation, we find 
that the helicity levels are filled according to 

. Ah=0) . (|/^| = 1) . (/^=0) . (|/^| = 1) . (\h\=2) _ 1 _ 1 _ 1 _ 1 _ 3 _ 6 

XqO ■ Xqi ■ Xqi ■ Xq2 ■ aq2 • aq2 ~ 9 ' 6 ' 6 ' 18 ' 18 ' 18 ' 

A few points about this result should be noted. Firstly, if we compare the helicity 



populations in (|3.2|) with their unaligned counterparts 

^ . Ah=G) . (\h\=i) . {h=o) . (|/i|=i) . A\h\=2) _ 1 . 1 . 2 _ 1 _ 2 _ 2 

aqo ■ Xqi ■ Xqi ■ Xq2 ■ X.Q2 ■ X.Q2 ~9'9'9'9'9'9' ^ ' 

we see that gluon fragmentation at high transverse momenta yields a sizable Xqj alignment. 
Secondly, the ratios in eqn. ( p.2| ) are consistent with the leading color-octet terms in the 



polarized g -^ Xa fragmentation functions which were calculated in ref. |T^ . The effect of 
the Xcj alignment upon the angular distribution of photons that result from the electric 
dipole transition Xa -^ J/il' + 7 as well as the induced J/^{j alignment were discussed in 



In order to completely determine the 0((as(2M(5)/7r)'^) contribution to ^f ^ Xqj fragmenta- 
tion, the interference between the graph in fig. |l|b and its one-loop corrections must be calculated 
along with the diagrams in fig. |l|a. Such interference terms have not been included so far in 
g -^ Xqj fragmentation function computations. 
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this previous article. We will not reiterate those findings of ref. [^ here, but we note that 
they may be simply recovered from spin symmetry arguments. Finally, uncertainty in the 
numerical value for the nonperturbative matrix element Hg drops out from the Xqj ratios 



in (3^). Therefore, they represent more reliable predictions than corresponding estimates 
for absolute g -^ Xqj fragmentation probabilities which depend upon the presently poorly 
constrained H^ parameter. 

We consider next gluon fragmentation to D-wave quarkonia. In the absence of a well- 
established nomenclature convention, we adopt the symbol Sqj to refer to L = 2 quarkonia 
with total angular momentum J. Quark model calculations indicate that n = 1 and n = 2 
dhj bottomonia lie below the BB threshold. On the other hand, all n = 1 Scj charmonia 
are predicted to have masses greater than 2Mo but less than M^ + Md* [ff9|-BT|. □ The 



l^Di and l^Ds states may decay to DD and are broad. But parity forbids their 1^D2 
and 1^D2 counterparts from decaying to two spinless mesons. The J = 2 charmonia are 
consequently quite narrow. Gluon fragmentation to 1^D2 paracharmonia has been found 
to yield a nonnegligible pp -^ 1^D2 + X cross section at the Tevatron [E^. Fragmentation 



to experimentally more accessible 1^D2 orthocharmonia starts at the same order in the 
velocity expansion but one lower order in as(2Mc)/7r. The 1^D2 states should therefore 
be produced in greater abundance. 

Gluon fragmentation to 5qj orthoquarkonia proceeds at 0{v'^) via Feynman diagrams 
like those illustrated in figs. 2a, 2b and 2c. The states emerging from the shaded ovals 
in these graphs correspond to the QQ pairs in the following Fock components of the 5qj 
wavefunction: 

I 5^, ) = 0(1) I QQ [^D'-p] ) + 0{v) I QQ [^Pf] g) + 0{v^) \ QQ p^f )] gg) + ---. (3.4) 

The diagrams in fig. [^a and 2b are suppressed relative to that in fig. ^jc by (Q;s(2MQ)/7r)^ 
and as{2MQ)/'K respectively. So the dominant contribution to g -^ Sqj fragmentation 
comes from the Fock component that contains the S-wave color-octet state. Its subsequent 
conversion to 6qj through long distance processes such as '^S[ -^ ^Dj + g + g and 
^5'} + g + g ^ ^Dj ' conserves jh = 1- We can thus determine the alignment which the 
5qj orthoquarkonium inherits from its transverse ^Sl progenitor by setting j^ = 0, j[ = 2 



Evidence for a 3.836 GeV 1^D2 charmoniuin state has recently been reported in ref. |11]. 



and L = 2 in eqn. ( |2.2|) . We find that the individual 5qJ hehcity states are populated 
according to 

.{h=0) . .{\h\ = l) . .(/^=0) . A\h\ = l) . .{\h\=2) . .(/^=0) . e(|/^| = l) . Mh\=2) . .(|/i|=3) 
"q1 • "q1 • "q2 • "q2 • "q2 ■ "q3 • "q3 • "q3 • "q3 

_373523710 15 (3-5) 

~50'50'30'30'30'75'75'75'75' 

Once 5qj orthoquarkonia are formed, spin symmetry fixes the leading order transfer 
of alignment from these D-wave states to S and P-wave orthoquarkonia descendants in 
subsequent strong and electromagnetic decays. Consider for example the decay mode 
5c2 -^ J 1"^ TTTT in the cc sector. This reaction conserves heavy spin angular momentum 
since the outgoing pions are soft. So we can again use eqn. (|2]^) to compute the alignment 
which the final J/t/j inherits from the initial 6c2'- 



J/^('^=o) = Br{6c2 -^ J/ip tttt) 



nr(/i=0) , 1 A\h\ = l) 2 (\h\=2) 

^"c2 + g "c2 + o "c2 

1 a('^=o) ^ ^ xd'^Ni) , 1 A\h\=2) 

-■-"02 "•" g "c2 "•" o "c2 



(3.6) 



The names of the helicity states in this expression denote fragmentation probabilities. 
Inserting the J = 2 D-wave results from eqn. (|3?5|), we deduce that dc2 -^ J/i/jtttt feeddown 
produces longitudinal and transverse J/tp^s in the ratio J/tp^^^^' : J/tp^'^'^^^ = 1 : 3.6. 

We now come to the most interesting application of spin symmetry ideas to heavy 
quarkonia systems. Specifically, we consider its implications for S-wave quarkonia for- 
mation. Prompt J/ip and ip' production at the Tevatron are currently the subject of 
significant theoretical and experimental interest. Previous studies have found that parton 
fusion alone underestimates the rate for J/ip production at large transverse momenta by 
more than a factor of 10. But if fragmentation is included, the disagreement between the- 



ory and experiment reduces to roughly a factor of 2 [|T5HT7|. The remaining discrepancy 
may be attributed to a number of theoretical uncertainties in perturbative QCD fragmen- 
tation function computations. The application of fragmentation ideas to J/tp production 
thus represents a qualitative success. 

The situation for prompt ip' production is totally different. The prediction for da{pp — > 
ip' + X)/dp± at high p± underestimates the measured cross section by a factor of 30 even 
after fragmentation is taken into account [^ . Recently, gluon fragmentation to some of the 
n = 2 Xcj states that lie above the DD threshold has been proposed as a possible significant 






source of ip'^s [18,24|. This mechanism requires an uncomfortably large radiative branching 



fraction Br{x'cj ~^ i^' +l) — 5-10% in order to explain the factor of 30 discrepancy between 
theory and data. However, a search for the photons that result from the electromagnetic 
transition along with an analysis of their angular distribution would at least provide a 
means for testing this proposal. 

There is another possible resolution to the t/j' surplus problem which does not in- 
volve any charmonia states above the DD threshold |2^. Instead, the answer may lie in 
fragmentation to higher Fock state components within the t/j' wavefunction 

\.l;') = 0{l)\cc['s['^]) + 0{v)\cc['Pf^]9) 

+ 0{v^) I cc pl^i^'^)] gg) + 0{v^) \ cc ['S^^^] g) + 0{v^) \ cc p^^^^] gg) + ■ ■ ■ . 

Gluon fragmentation to the ccf^S"} ] pair in the leading Fock component starts at 
0(f'^(as(2Mc)/7r)'^) via Feynman diagrams like the one illustrated in fig. ^. Fragmen- 
tation to the cc states in all other Fock components of \i^') takes place at higher order 
in the velocity expansion. But the heavy charm-anticharm pairs in the subleading Fock 
components can be formed at lower orders in as(2Mc)/7r. In particular, the transforma- 
tion shown in fig. |]c of a hard gluon into a ^S"} cc pair followed by its conversion into 
-(/;' via the emission or absorption of two soft gluons occurs at 0(f^as(2Mc)/7r). If we 
evaluate the relative strengths of the -i/;' production processes in figs. 2a and 2c by set- 
ting v"^ ~ 0.3 and as(2Mc)/7r ~ 9.0 x lO'^, we find v^{as{2Mc)/T:f ^ 1.2 x 10"^ and 
f ^(as(2Mc)/7r) ^ 1.4 x 10""^. Hence, gluon fragmentation to the color-octet cc pair within 
the last Fock state of eqn. (|3.7| ) may naively be an order of magnitude more important 
than that to the color-singlet pair within the first Fock state. If so, its inclusion could 
resolve the ip' problem. Q 

Without knowing the value for the nonperturbative matrix element associated with 
g -^ cc^S]^ ] fragmentation, we cannot draw any definite conclusion regarding this possible 
'0' production mechanism. But spin symmetry provides a powerful means for testing 
the hypothesis that this proposal accounts for most of the prompt ■(/'"s observed at the 
Tevatron. Recall that the "^Sl cc pair in fig. ||c is transverse up to 0{M^,/qQ) corrections. 
This intermediate ji = state transforms at leading order in the velocity expansion into a 
j'^ = t/j' through the L = emission or absorption of two soft gluons. The outgoing i/;' in 



^ Gluon fragmentation to the cc pair within the |cc[^S'J ]gg) Fock component of jV') also 
enhances prompt J/tjj production. However, its impact is smaller than that for ip' production 
since high transverse momenta J/ip^s predominantly come from g -^ Xa fragmentation. 
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fig. I^c has exactly the same angular momentum quantum numbers as the incoming hard 
gluon. Therefore, tj^'^s produced at high transverse momenta via gluon fragmentation to 
the final Fock state in eqn. ( |3.7| ) are 100% transversely aligned! 

This striking prediction can be tested by measuring the lepton angular distribution in 
the if:' —>■ i~^£~ decay channel. The decay products are generally distributed according to 

"^^ ii;' ^i+r)oc{l + acos^O) (3.8) 



dcos9 



where 9 denotes the angle between the lepton three-momentum in the t/^' rest frame and the 
t/j' three-momentum in the lab frame. The parameter a ranges over the interval — 1 < a < 1 
and equals unity for totally transverse i/^'^s. Subleading corrections to the infinite mass 
limit prediction will shift the value for a to less than one. But experimental observation 
of a lepton angular distribution close to dV/d cos 9 oc (1 + cos^ ^) would provide strong 
support for this resolution to the '0' surplus problem. 
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Figure Captions 

Fig. 1. Feynman diagrams wiiidi mediate gluon fragmentation to P-wave Xqj ortho- 
quarkonia at order (a) t'^(as(2MQ)/7r)^ and (b) v^as{'2.MQ)/'K. 

Fig. 2. Feynman diagrams which mediate gluon fragmentation to D-wave 5qj ortho- 
quarkonia at order (a) v'^ {as{2MQ)/Tif , (b) v'^ {as{2MQ)/TxY and (c) v'^ as{2MQ)/TX. 
These same diagrams also mediate gluon fragmentation to S-wave iI)q states at 
order (a) v'^{as{2MQ)/Tif, (b) v'^{as{2MQ)/TiY and (c) i'^as(2Mg)/7r. Crossed 
graphs in (a) and (b) are not displayed. 
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